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ABSTRACT 

Validity  of  seven  analytical  formula  as  well  as  linear  and  non-linear  mul¬ 
tiple  regressive  schemes  was  tested  using  field  data  from  the  Outer  Banks,  North 
Carolina.  Generally,  agreement  proved  unoatiafactory.  Field  experiences  in¬ 
dicate  that  the  longshore  current  is  a  velocity  field  consisting  of  a  multitude  of 
velocity  vectors  whose  basic  pattern  varies  depending  upon  the  regimes  of  wave- 
current- topography  interaction.  The  need  to  recognize  topography  as  a  respond¬ 
ing  variable  as  well  as  a  process  variable  in  the  physical  scheme  of  longshore 
current  is  emphasized.  ^ 


INTRODUCTION 

One  may  quote  today  more  than  a  dozen  longshore  current  formulae  which 
have  been  developed  on  the  basis  of  simplified  analytical  schemes.  Table  1  lists 
the  exponents  of  individual  variables  when  these  formulae  are  reduced  to  a  pro¬ 
duct  form.  Scanning  down  each  column,  it  la  found  that  these  formulae  contain 
widely  vaired  contributions  of  the  component  variables.  For  instance,  note  the 
incidence  angle,  ©.  This  Important  variable  is  taken  aea  sine  function  in  some 
formulae,  as  a  cosine  function  in  others,  with  0  or  20,  carrying  different  ex¬ 
ponents.  Thus,  the  question  arises  as  to  which  of  these  formulae  gives  the  best 
approximation  to  phenomena  as  observed  In  the  field.  Also.  It  has  not  been  estab¬ 
lished  whether  the  simplifying  assumptions  used  in  deriving  these  forumlae  would 
give  a  marginal  or  a  significant  distortion  of  reality.  These  problems  are  dealt 
with  in  this  paper  on  the  basis  of  field  experience. 


Three  different  procedures  of  quantitative  analysis  are  applied.  The  first 
attempts  to  approximate  the  mean  alongshore  components  of  longshore  current 
velocities  through  a  linear  combination  of  individual  variables;  a  method  known 
as  multiple  linear  regression.  The  second  is  a  non- linear  approximation  using  a 
product  form  of  optimum  powers  of  individual  variables.  By  a  simple  mathemati¬ 
cal  manipulation  this  form  la  reduced  to  the  case  of  linear  regression.  ^Thf  th^rdj 
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tests  the  goodness  of  fit  between  formula  -  derived  predictions  and  actual  field 
observations.  The  formulae  tested  include  those  by  Putnam -Munk- Traylor 
(1949),  Inman-Quinn  (1951),  Nagal  (1954),  Brebner-  Kamphu  s  (1963),  Galvin- 
Eagleson  (1965),  and  Inman- Bagnold  (1962). 

The  second  half  of  this  paper  deals  with  realistic  recognition  of  the  role  of 
topography  in  the  physical  scheme  of  longshore  currents.  The  longshore  current 
is  recognized  as  a  velocity  field  which  consists  of  a  multitude  of  velocity  vectors. 
Field  experiences  are  presented  indicating  that  the  topography  is  a  responding 
variable  as  well  as  a  process  variable,  and  that  depending  upon  the  regimes  of 
wave- current -topography  interaction  some  meaningful  patterns  of  velocity  field 
could  take  place. 


DATA  BACKGROUND 

The  data  utilized  in  this  paper  include:  (I)  a  series  of  field  measurements 
conducted  by  the  Coastal  Studies  Institute,  Louisiana  State  University,  on  the 
Outer  Banks  beach,  North  Carolina:  (2)  field  and  laboratory  data  from  Putnam 
Munk-Traylor;  and  (3)  field  and  laboratory  data  from  scattered  sources. 

The  field  data  an  the  Outer  Banks  beach  (Figure  1)  were  obtained  by  meas¬ 
uring  the  current  velocities  four  times  within  each  tidal  cycle,  simultaneously 
with  other  related  variables  such  as  winds,  waves,  tides,  air  and  water  tempera¬ 
tures,  beach  profiles  and  sediment  characteristics.  Since  the  coordinated  cover¬ 
age  of  all  these  variables  required  the  use  of  a  stable  platform,  the  field  activi¬ 
ties  were  maintained  close  to  a  fishing  pier  (Figure  2).  The  current  velocity  was 
measured  by  timing  the  movement  of  a  dye  patch  over  a  fixed  distance  of  70  feet 
in  either  direction  from  this  pier.  The  wave  data  were  supplied  by  the  Btep- 
resistance  gage  of  the  Coastal  Engineering  Research  Center,  which  was  main¬ 
tained  alongside  the  pier  at  a  position  approximately  5  meters  (15  feet)  deep.  The 
field  operation  was  continued  for  approximately  6  months  between  December, 

1963  and  May,  1964. 

QUANTITATIVE  TESTS  OF  SIMPLIFIED  SCHEMES 
MULTIPLE  LINEAR  REGRESSIVE  SCHEME 

The  analytical  technique  of  a  multiple  linear  regressive  scheme  has  been 
excellently  described  by  Harrison  and  Krumbein  (1964).  The  basic  equation  is 
written: 


Y  =  bQ  +  bj  Xi  +  t>2  X2  +  .  .  .  +  bfc  Xk, 

in  which  Y  is  the  dependent  variable,  bQ,  bj ,  .  .  .  the  partial  regression  coef¬ 
ficients,  and  Xj,  Xg,  .  .  .  the  independent  variables.  In  the  present  case,  the 


Fig.  1.  Location  of  study  area;  Nags  Head,  the  Outer  Fig.  -.  Field  instrumentation  and  beach  profile  at  study  site. 
Banks,  North  Carolina. 
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dependent  variable  is  the  mean  alongshore  component  of  the  longshore  current 
velocities.  The  independent  variables  are  chosen  as  follows:  two  measures  of 
wave  height  (H,  wave  height  measured  by  the  CERC  wave  gage;  and  H^,  breaker 
height  visually  observed  over  the  inner  bar),  the  wave  period  (T),  the  angle  of 
incidence  (6  or  sin  26,  the  angle  between  the  wave  creBt  and  the  shoreline),  the 
mean  slope  of  the  surf- zone  bed  (M),  and  the  alongshore  component  of  wind  velo¬ 
cities  (W).  The  signs  are  positive  for  V,  6  and  W  arriving  from  the  north,  and 
negative  for  these  arriving  from  the  south. 

In  the  numerical  computation,  a  special  IBM  program  was  designed,  ao 
that  all  the  possible  combinations  of  different  numbers  of  independent  variables 
could  be  dealt  with.  The  relative  contribution  of  each  independent  variable  to 
observed  mean  longshore  current  velocities  was  evaluated  by  noting  the  reduction 
in  R2  between  two  regressions  including  and  excluding  this  particular  variable. 

In  statistical  terminology,  R2  Is  a  measure  of  the  fraction  of  the  total  variance  of 
the  dependent  variable  which  is  accounted  for  by  the  regression.  The  square  root 
of  this  measure,  R,  is  the  multiple  correlation  coefficient.  The  results  of  the 
multiple  linear  regression  analysis  are  shown  in  Table  2. 

The  analysis  shows  that  a  considerably  high  level  of  explained  variation 
(up  to  72  per  cent  in  R2)  could  be  attained  by  use  of  multiple  linear  regression. 
The  most  Important  single  variable  affecting  the  mean  longshore  current  velocity 
turns  out  to  be  the  angle  of  wave  Incidence,  6  or  sin  26,  which  alouc  accounts  for 
approximately  <J8  per  cent  out  of  the  maximum  R2  level  of  72  per  cent.  This  is 
followed  by  ihe  wind  velocity  (accounting  for  the  mere  fraction  of  3. 1  per  cent  in 
R2),  the  wave  height  (I.  3  per  cent  in  R2),  the  bed  slope  (0. 7  per  cent  in  R2),  and 
the  wave  period  (0.0  per  cent  in  R2)  in  the  order  named.  It  appears,  therefore, 
that  so  far  as  the  Outer  Banks  data  ere  concerned,  the  mean  longshore  current 
velocities  can  be  explained  to  a  reasonably  high  degree  through  a  simple  linear 
regression  containing  the  angle  of  incidence  as  the  sole  Independent  variable,  as 
follows: 


V  -  -0.264  +  2.958  sin  26, 

6  in  degrees,  and 
V  in  ft.  / sec. 

Similar  results  were  obtained  by  using  the  field  and  laboratory  data  of 
Putnam- Munk-Traylor  (1949)  (Table  3).  The  levels  of  explanation,  R2,  using 
these  data  were  89.1  and  65.3  percent  for  the  laboratory  and  field  obser¬ 
vations,  respectively.  In  both  of  these  sets  of  data,  the  most  important  single 
variable  was  also  the  angle  of  incidence,  which  accounted  for  81.5  per  cent  of 
the  total  R2  of  89. 1  per  cent  and  47.  9  per  cent  out  of  the  total  R2  of  65.  3  per 
cent,  respectively. 

These  results  are  conflicting  with  those  previously  reported  by  Harrison 
and  Krumbein  (1964),  who  used  data  obtained  on  the  Virginia  beach,  approximately 
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70  miles  north  of  our  study  area.  The  angle  of  wave  incidence,  which  proved  to 
be  the  most  influential  variable  both  in  the  Outer  Banka  and  the  Putnam- Munk- 
Yraylor  data,  ranked  an  the  fifth  most  important  variable  in  their  analysis.  The 
most  important  variable  turned  out  in  their  analysis  to  be  the  wave  period,  which 
proved  to  be  insignificant  in  our  analysis. 

The  discrepancies  just  described  may  probably  indicate  the  very  limim- 
tions  irherent  to  the  multiple  linear  regressive  scheme.  This  scheme  only  pro¬ 
vides  a  linearized  description  of  a  certain  physical  scheme  through  a  particular 
set  of  sample  data  available  for  analysis.  Consequently,  depending  upon  the  part¬ 
icular  samples  analyzed  as  well  as  the  non-linearity  characteristics  of  the  ori¬ 
ginal  physical  scheme,  a  considerable  distortion  of  reality  might  result.  In  the 
present  cape,  the  sample  characteristics  may  relate  to  the  regional  and  seasonal 
regimes  of  the  beach-ocean-atmosphere  interaction  system.  While  both  the  Outer 
Banks  data  and  the  Putnam- Munk- Traylor  data  was  obtained  during  single  seasons, 
the  former  during  January  to  March  and  the  latter  during  February  to  April,  those 
of  Harrison-Krumbein  were  taken  from  scattered  periods  between  February  and 
July. 

MULTIPLE  NON-  LINEAR  REGRESSIVE  SCHEME 

Review  01'  various  published  formulae  suggests  that  the  relationship  be¬ 
tween  the  mean  longshore  current  velocity  and  the  independent  variables  might  be 
approximated  by  a  non-linear  relationship,  such  as  a  product  form  of  Independent 
power  variables,  namely 

V  =  e^X^X^2  .  .  .  x£k 

By  taking  the  logarithm  on  both  sides,  one  obtains 


log  V  =  fcQ  +  bx  log  Xj  +  b2  log  Xjj  +  .  .  .  +  bfc  log  Xk. 


By  regarding  the  logarithms  of  individual  variables  as  the  elementary  variables 
the  problem  reduces  to  that  of  a  multiple  linear  regressive  scheme.  The  results 
of  analysis  are  summarized  in  Table  4. 

In  the  Outer  Banks  data,  the  multiple  correlation  level  based  on  this  non¬ 
linear  scheme  is  found  to  be  appreciably  low,  i.  e.  51  per  cent.  Also  unlike  the 
results  of  linear  regression,  the  most  important  variable  affecting  the  mean  long¬ 
shore  current  velocity  proves  to  be  the  wave  height,  accounting  for  approximately 
22  per  cent  out  of  the  total  25  per  cent  in  R2.  The  angle  of  incidence  is  the  least 
important  of  all  the  independent  variables.  The  F-tests  show  that  any  combination 
of  variables  excluding  either  the  breaker  height  or  the  shallow-water  wave  height 
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(measured  by  the  CERC  gage)  was  statistically  insignificant  at  the  1  per  cent  level 
of  probability. 

A  similar  analysis  using  the  laboratory  and  field  data  of  Putnam- Munk- 
Traylor  indicates  a  comparable  level  of  correlation,  approximately  61  per  cent 
and  82  per  cent,  respectively  (Table  5).  Interestingly,  the  importance  of  wave 
height  i6  noted  only  In  the  laboratory  data  (43  per  cent  in  ft2).  The  angle  of  inci¬ 
dence,  which  was  insignificant  in  the  non-linear  regression  based  on  the  Outer 
Banks  data,  proves  to  be  important  in  both  sets  of  the  Putnam-Munk-Traylor 
data,  accounting  for  24  per  cent  and  45  per  cent  in  R2,  respectively. 

Consequently,  it  is  stated  that  the  non-linear  regression  of  a  product  form 
does  not  necessarily  improve  the  level  of  correlation.  Comparing  between  the 
linear  and  non-linear  regressive  schemes  and  also  between  the  data  from  th<‘  dif¬ 
ferent  sources  analyzed  here,  a  considerable  confusion  arises  as  to  the  relative 
contribution  of  individual  variables.  Apparently,  further  study  is  needed  before 
the  multiple  regressive  schemes  become  an  effective  means  of  prediction. 

TEST  OF  FORMULAE 

General  -  All  the  known  formulae  (Table  1)  have  been  derived  from  simpli¬ 
fied  physical  schemes  in  which  the  bottom  topography  is  replaced  by  an  inclined 
plane.  The  approach  by  Inman  and  Bagnold  (1962),  Shadrin  (1961)  and  Bruun 
(1963)  has  further  assumed  the  presence  of  a  longshore  bar  or  multiple  bars  with 
or  without  regular  gaps  along  the  bar  crests.  Theoretical  derivation  of  physical 
relationships  assumes  that  a  steady  state  of  equilibrium  is  established  in  terms  of 
either  momentum  exchange,  energy  conservation,  mass  continuity  or  surface 
gradient.  The  mathematical  expressions  of  the  known  analytical  formulae  are 
summarized  in  Table  1. 

Putnam-Munk-Traylor  formulae  (1949),  revised  by  Inman-Quinn  (1951)  - 
Putnam-Munk-Traylor  were  the  first  to  introduce  a  simplified  physical  scheme 
amenable  for  theoretical  analysis.  Using  the  extensive  field  observations  on  the 
California  coasts,  Inman  and  Quinn  later  demonstrated  that  the  friction  coefficient, 
K,  contained  in  the  original  Putnam-Munk-Traylor  formulae  could  be  revised  in 
accordance  with  the  relationship: 


K  oC  V1-5 


Computation  using  the  Outer  Banks  data  indicates  tha*.  the  prediction  based  on  the 
revised  momentum  formula  exceeds  the  observed  current  velocities  by  the  factor 
of  2.  67  (Figure  3).  The  friction  coefficients  computed  from  the  original  Putnam- 
Munk-Traylor  formula  (Figure  4)  further  indicate  that  the  friction  coefficient  may 
not  be  as  definite  a  function  of  the  current  velocities  as  suggested  by  Inman  and 
Quinn.  An  adequate  explanation  of  this  difference  between  the  results  of  Inman 
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Predicted  Velocity,  Ft/Sec 

Fig.  3.  Prediction  using  Inman-Quinn- revised  Putnam-Munk-Traylor  momentum 
formula  versus  field  observations  from  (X>ter  Banks,  N.C. 


Fig.  4.  Friction  coefficients  computed  from  Putnam-Munk-Traylor  formula  and 
Nagai  formula  versus  current  velocities. 
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a  rid  Quinn  and  the  present  authors  in  not  immediately  possible.  However,  it  may 
be  that  the  K  value  is  influenced  not  only  by  the  bed  roughness  but  also  by  the 
effects  of  other  energy  dissipative  sources  •  namely  the  internal  turbulence  in  the 
surf  zone,  and  the  sediment  movement  and  the  topographic  changes  which  take  up 
a  considerable  portion  of  wave  energy  contributed  here. 

Naoai  formula  (1954)  -  Referring  to  the  study  of  Housely  and  Taylor  (195?), 
Shadrin  (1951)  questioned  the  application  of  the  solitary  wave  theory  to  approxi¬ 
mate  breaker  characteristics,  such  as  was  used  in  the  derivation  of  the  Putn&m- 
Munk- Traylor  formulae.  This  problem  can  be  checked  through  a  formula  pro¬ 
posed  by  Nagai  (19  5  4),  whose  derivation  used  the  Airy's  first-order  approxima¬ 
tion  along  with  the  same  physical  scheme  as  proposed  by  Putnam-Munk- Traylor 
(Table  1). 

Figure  4  includes  the  friction  coefficients  computed  from  this  formula 
using  the  Outer  Banks  data.  The  scatter  of  the  plots  suggests  that  the  application 
of  the  oscillatory  wave  theory  does  not  substantially  improve  the  goodness  of  fit 
with  observation. 

Brebner-Kamphulw  formula  (1963)  -  The  energy  and  momentum  formula 
proposed  by  these  Investigators  contain  the  deep-water  equivalents  of  wave  char¬ 
acteristics  on  the  consideration  that  these  are  die  fundamental  variables  available 
through  wave  fores  sting  techniques. 

The  prediction  using  the  energy  and  momentum  formulae  exceede  the 
observed  velocities  by  the  factor  of  1. 76  and  2. 73,  respectively  (Figure  5).  The 
disagreement  is  too  obvious  to  be  accounted  for  by  the  fact  that  our  observed  velo¬ 
cities  tended  t  >  be  small  due  to  the  proximity  of  a  fishing  pier. 

Galvin- Eagleson  formula  (1965)  -  Figure  6  shows  the  South  flowing  and 
the  Nr  ifh  flowing  currents  separately.  With  both  of  these  currents  combined,  the 
prediction  exceeds  the  field  observation  by  the  factor  of  3.  42. 

Inman- Bagno  Id  formula  (1962)  -  The  basic  scheme  of  derivation  assumes 
that  the  mass  of  water  released  by  the  breaker  is  preserved  within  the  surf  zone 
until  it  is  channeled  out  through  the  rip  outflow  at  a  downstream  position.  Thus, 
the  formula  contains  two  unknowns:  the  fraction  of  the  total  wave  mass  which  is 
contributed  to  the  current,  denoted  by  _s,  and  the  Interval  between  adjacent  rips, 
denoted  by  _l  (Table  1).  Using  the  Outer  Banks  data,  the  terms  excluding  these 
unknowns  are  computed  and  plotted  against  observed  velocities  in  Figure  7.  The 
plots  suggest  the  relationship: 


S'  1  :  30  (  unit  in  feet  ) 


Assuming  that  the  s  value  ranges  between  2  and  10  per  cent  (Galvin- Eagleson, 
1965;  Brebaer-Kamphuls,  1963;  Inman- Bagnold,  1962),  the  rip  interval,  1,  should 
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Predicted  Velocity.,  P^ec 

Fig.  5.  Prediction  using  Brebner-Kamphuis  formulae  versus  field  observations 
from  Outer  Banks,  North  Carolina. 


Fig.  6.  Prediction  using  Galvin-Eagleson  formula  versus  field  observations 
from  the  Outer  Banks,  N.C. 
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fall  between  1050  ft.  to  210  ft.  -  the  order  roughly  comparable  to  the  case  of  the 
Outer  Banks  beach  as  well  as  that  previously  reported  by  Shepard  (1950). 


VARIABILITY  OF  LONGSHORE  CURRENTS 

GENERAL 

It  must  be  mentioned  that  the  tests  thus  far  described  utilized  selected  data, 
excluding  anomalies  such  as  the  currents  opposed  to  the  direction  of  wave  ap¬ 
proach,  the  non- zero  current  velocities  with  perpendicular  wave  incidences,  the 
zero  current  velocities  with  oblique  incidences,  and  the  wind  directions  opposed  to 
the  current  directions. 

In  all,  these  anomalies  occurred  in  35  per  cent  of  the  field  observations, 
implying  that  the  formulae  based  on  simplified  schemes  failed  to  explain  a  consid¬ 
erable  portion  of  the  actual  phenomenon.  Even  with  the  selected  data,  the  dis¬ 
agreement  between  the  predicted  and  the  observed  velocities  was  appreciable.  It 
may  be  that  while  the  formulae  are  designed  to  give  the  mean  velocity  in  the  surf 
zone,  no  known  method  of  observation  is  capable  of  yielding  an  unbiased  mean,  or 
that  the  simplified  analytical  schemes  fail  to  account  for  the  actual  physical  mec¬ 
hanism  of  the  longshore  current  phenomenon. 

From  the  hydrautics  point  of  view,  the  longshore  current  is  a  case  of  an 
unsteady  flow  confined  in  a  time-variant  boundary.  The  current  field  contains 
velocity  vectors  which  are  variable  both  in  time  and  space,  and  the  bottom  topo¬ 
graphy  is  changing  constantly.  Generally  speaking,  the  temporal  variability  may 
be  attributed  to  the  stochastic  nature  of  the  input  waves,  and  the  spatial  variabili¬ 
ty  to  the  interplay  between  waves,  currents  and  the  topography.  It  is  likely  that 
this  latter  effect  is  a  reciprocal  process,  in  which  a  change  of  each  variable  is 
fed  back  to  the  behaviors  of  other  associated  variables.  Consequently,  although 
the  simplified  analytical  schemes  assume  the  topography  to  be  a  fixed  boundary 
and  the  waves  and  the  currents  to  be  steady  in  time,  the  actual  phenomena  must 
be  comprehended  in  view  of  the  dynamic  regimes  of  interaction  between  these  vari¬ 
ables. 

TEMPORAL  VARIABILITY 

Owing  to  the  difficulty  of  field  measurement,  quantitative  data  pertaining 
to  the  temporal  variability  are  particularly  scarce.  The  indication  of  the  tempo¬ 
ral  variability  of  longshore  currents  has  been  >«port«J  by  Shepard  and  Inman 
(1950)  on  the  Scripps  beach,  La  Jolla,  California.  Pulsation  of  the  rip  outflow 
observed  on  this  beach,  with  tne  average  recurrent  period  of  J.7  minutes,  has 
been  reported  to  coincide  with  that  of  the  sur?  beat  activities.  According  to  obser¬ 
vations  on  the  Niigata  beach,  Japan  (Fujiki,  1957),  the  rip  pulsation  was  appar¬ 
ently  associated  with  *.he  angle  of  wave  incidence;  the-  pulsation  gave  way  to  a  con¬ 
tinuous  diagonal  outflow  ao  the  obliquity  of  wave  incidence  increased. 
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Fig.  7.  Inman- Bagnold  formula  tested  by  checking  the  magnitude  of  mass  con¬ 
tribution  to  currents  by  waves,  s,  and  rip  separation,  l. 


Fig.  8.  Coefficient  of  variation  versus  mean  of  10-minute  average  velocities. 
Data  from  Niigata,  Japan  (1958,  a). 
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The  data  on  the  Niigata  beach,  Japan,  obtained  with  a  self-contained  auto¬ 
matic  recorder  (ONO-Meter),  probably  represent  the  largest  source  of  informa¬ 
tion  regarding  this  subject  (Niigata,  1958, a, b).  According  to  these  observations, 
the  temporal  variability  beyond  the  order  of  10-minute  period  appears  to  be  rela¬ 
tively  low.  The  standard  deviation  of  the  10-minute  average  velocities  remained 
below  0.  3  ft/ sec  throughout  the  entire  runs  of  data.  This  is  particularly  meaning¬ 
ful  when  one  considers  that  during  the  period  of  current  measurements  breakers 
as  high  as  15  ft  were  observed.  The  data  by  Inman  and  Quinn  (1951),  taken  during 
swell  activities  and  probably  containing  the  effects  of  both  temporal  and  spatial 
variabilities,  have  resulted  in  standard  deviation  of  0.5  ft/ sec.  In  the  Niigata 
data,  the  coefficient  of  variation  of  the  10-minute  average  velocities  was  merely 
0.  3  for  the  velocities  near  1. 0  ft/sec  (,  3m/sec),  decreasing  gradually  to  as  low 
as  0. 1  for  both  lower  and  higher  velocities  (Figure  6).  According  to  Inman  and 
Quinn,  the  coefficient  of  variation  remained  not  lower  than  1.0  for  all  the  veloci¬ 
ties  observed.  Indications  are  that  the  proportion  of  the  scatter  in  current  veloci¬ 
ties  which  is  attributable  to  the  temporal  variability  is  smaller  than  to  the  areal 
variability.  Interestingly,  current  reversal  was  reported  to  occur  twice  as  fre¬ 
quently  near  the  offshore  ter  as  inshore  of  this  position.  In  general,  the  tempo¬ 
ral  variability  appears  to  be  more  pronounced  near  the  offshore  bar  than  inaide 
the  surf  zone. 

SPATIAL  VARIABILITY 

The  spatial  variability  of  longshore  currents  perpendicular  to  the  shore¬ 
line  has  been  recognized  in  laboratory  models  (Galvin- Eagle  son,  1964;  Shiniano, 
et  al,  1957)  as  well  as  in  the  field  (Ajbulatov,  et  al,  1966).  Variability  in  the 
alongshore  direction  has  received  much  attention  in  recent  years  (Inman- Bagnold, 

1962) .  This  type  of  variability  includes  the  processes  of  gradual  acceleration,  de¬ 
celeration,  stagnation,  and  reversal,  and  has  even  been  recognized  in  a  labora¬ 
tory  model  with  straight  contours  (Galvin- Eagieson,  1964;  Brebner-Kamphuis, 

1963) .  Variability  with  depth  is  probably  the  least  known  feature  but  highly  im¬ 
portant  in  view  of  Itc  implications  on  sediment  transport.  Actual  observation 

has  been  made  mostly  bey  .nd  the  longshore  bar  (Shepard -Inman,  1950;  Sonu,  1961), 

Certain  systematic  combinations  of  these  variability  features  seem  to 
emerge  depending  upon  the  interplay  among  the  associated  variables.  The  regu¬ 
larly  spaced  rip  currents  (McKenzie,  1958;  Shepard,  1950)  represent  one  of  such 
examples.  Accordingly,  the  longshore  current  velocity  field  may  be  classified 
into  four  distinctive  categories  in  the  light  of  the  regimes  of  this  interplay  -  the 
hypothetical  regime,  the  natural  equilibrium  regime,  the  transitional  regime,  and 
the  forced  equilibrium  regime. 

HYPOTHETICAL  REGIME 

The  topography  is  a  fixed  bed,  and  the  dynamic  Interaction  exists  only 
between  currents  and  waves.  Although  the  case  represented  by  this  regime  may 
seldom  arise  under  natural  conditions,  except  when  the  beach  is  very  steep  and 
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•  H-1.5M  T>8$ec  8  s  +5® 

o  H»2.0M  Tx  7 s«c  9--15* 


Fig.  9.  Longshore  current  velocity  field  associated  with  sand  waves.  Upper 
diagram  for  H“1.6m,  T"8  sec.  and  8“+50  (from  right  of  normal  to 
the  shore).  Lower  diagram  shows  mean  alongshore  velocity  compon¬ 
ents  at  200  m  from  the  shore,  as  indicated  by  thick  arrows  In  the 
upper  diagram.  (Adapted  from  Inokuchi,  i960). 


Fig.  10.  Currents  stagnate  over  the  shoal,  describing  circulating  paths.  Data 
■  from  Niigata,  Japan  (1958, b). 
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the  wave  action  very  weak,  the  hydraulic  mechanism  revealed  in  this  simplified 
case  appears  to  render  information  of  basic  interest. 

Laboratory  experiments  using  a  fixed  slope  (Galvin- Eagle  son,  1965)  and  a 
sand  bed  (Shimano  et  al,  1957}  have  demonstrated  the  presence  of  a  velocity  distri¬ 
bution  across  the  surf  zone,  with  the  maximum  velocity  occurlng  at  positions  ap¬ 
proximately  10  to  40  per  cent  of  the  surf-zone  width  from  the  shoreline.  It  has 
also  been  reported  (Galvin- Eagleson,  1965)  that  "the  width  of  the  longshore  cur¬ 
rent  expands  in  the  direction  of  flow  downstream  of  ti„j  obstacle."  ThiB  observa¬ 
tion  is  consistent  with  a  prediction  resulting  from  the  continuity  approach  (Inman- 
Bagnold,  1962)  that  the  distance  between  the  bar  crest  and  the  shoreline  must  in¬ 
crease  in  the  direction  downstream  from  the  point  of  a  rip  current. 

NATURAL  EQUILIBRIUM  REGIME 

The  natural  equilibrium  regime  represents  the  case  in  which  the  inter¬ 
action  between  waves,  currents,  and  the  topography  has  reached  a  state  of  dynam¬ 
ic  equilibrium.  Unlike  in  the  hypothetical  regime,  the  topography  participates  in 
this  interaction  as  a  responding  variable  as  well  as  a  process  variable.  Of  parti¬ 
cular  interest  is  the  fact  that,  when  acted  upon  by  longshore  currents  over  a  sus¬ 
tained  period  of  time,  the  nearshore  bed  develops  a  series  of  sand  waves  having 
elongated  ridges  and  troughs  directed  at  angles  to  the  shoreline  (Sonu  and  Russell, 
1966).  These  sand  waves  are  subsequently  reworked  by  waves  and  transformed 
into  an  alternate  sequence  of  lunate  bars  and  intermediate  shoals,  with  a  shoal- 
to- shoal  separation  of  several  hundred  to  several  thousand  feet  alongshore.  The 
longshore  current  velocity  field  corresponding  to  this  type  of  rhythmic  topography 
(Hom-ma  and  Sonu,  1962)  takes  on  a  typical  pattern  of  natural  equilibrium  re¬ 
gime.  Figure  9  shows  the  laboratory  reproduction  of  such  a  velocity  field  mod¬ 
eled  after  the  Niigata  beach,  Japan  (Inokuchi,  1960).  Although  the  current  pat¬ 
tern  is  a  group  of  circulating  cells  associated  with  the  shoals  and  lunate  bars, 
the  alongshore  velocity  components  plotted  along  the  shoreline  give  only  alternate 
sequences  of  converging  and  diverging  currents.  Figure  10  represents  the  actual 
field  observation  of  circulating  water  movements  in  the  vicinity  of  a  shoal,  on  the 
Niigata  beach,  Japan.  The  alongshore  velocity  components  are  generally  greater 
in  front  of  the  embayment  than  near  the  apex  or  the  shoal,  the  latter  being  associ¬ 
ated  frequently  with  stagnation  or  even  reversal  of  current  directions.  Figure  11 
shows  the  breaker  distribution  over  the  rhythmic  topography,  taken  from  the 
Niigata  beach.  The  breaker  type  is  diversified  into  the  spilling  type  over  the 
shoal  and  the  plunging  typo  over  the  bar. 

Figure  12  shows  the  longshore  current  velocity  field  actually  observed  by 
the  present  investigators  on  the  Outer  Banks  beach.  Waves  with  a  significant 
height  of  2.  5  feet  and  a  period  of  7. 8  seconds  arrived  at  30*  to  the  shore,  and 
there  was  no  wind.  The  waves  broke  by  spilling  over  the  shoal,  ri  leasing  a  con¬ 
siderable  amount  of  water  mass  directed  onshore.  This  free  mass  of  water  was 
further  joined  by  the  backwash  from  the  subaerial  beach  and  gradually  displaced 
downstream  by  the  alongshore  component  of  oblique  wave  drags.  Reaching  the 
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Fig.  11.  Breakers  over  lunate  bars  and  shoals,  characterizied  by  plunging  and 
spilling  types,  respectively.  Photo  showing  the  Niigata  beach,  Japan, 
taken  by  Kokusal  Aerial,  Inc.,  Tokyo,  October  12,  1955. 


Fig.  12.  Longshore  current  velocity  field  observed  on  the  CXiter  Banks,  North 
Carolina,  October  24,  1965. 
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center  of  the  jmbayment,  the  water  movement  became  stagnant,  due  apparently  to 
the  shoreline  curvature  which  resulted  in  a  decrease  in  the  alongshore  component 
of  wave  drag.  The  continuous  influx  of  water  caused  a  hydrostatic  potential  to 
build  up  at  this  position,  which  was  eventually  released  toward  offshore  as  the 
mass  outflow  of  a  rip  current.  This  phenomenon  is  predominantly  a  process  of 
mass  transfer.  On  the  other  hand,  the  waves  broke  over  the  bar  crest  by  plunging, 
but  much  of  their  energy  was  absorbed  into  secondary  waves  which  were  regener¬ 
ated  above  the  trough.  Consequently,  no  appreciable  amount  of  onshore  mass 
movement  occurred  beyond  the  trough,  the  water  released  by  the  breaker  under¬ 
going  a  gradual  displacement  downstream  along  a  relatively  constricted  path  over 
the  bar  crest.  Surprisingly,  the  dye  patch  entering  this  zone  created  a  clear  de¬ 
marcation  with  the  water  offshore  and  was  never  diffused  into  the  area  between 
the  bar  and  the  shoreline.  This  phenomenon  is  predominantly  a  process  of  momen¬ 
tum  transfer.  Thus,  it  appeared  that  the  current  field  just  described  consisted  of 
both  the  continuity  and  the  momentum  processes  operating  simultaneously.  Evans 
(1939)  has  reported  a  similar  case  from  the  Lake  Michigan  shore. 

The  longshore  current  velocity  field  associated  with  the  natural  equilibrium 
regime  may  be  a  more  general  occurrence  than  is  usually  believed.  The  rhyth¬ 
mic  topography  associated  with  this  regime  have  been  reported  from  various  sour¬ 
ces  of  the  world,  including  the  Azov  and  Caspian  Sea  coasts  (Shadrin,  1961); 
Egorov,  1951;  Kashechkin-Uglev),  the  French  Mediterranean  coast  (Riviere  et  al, 
1961),  the  Italian  coast  (King-Williams,  1949),  the  Danish  North  Sea  coast  (Bruun, 
1954),  the  Japanese  coast  (hom-ma-Sonu,  1962;  Mogi,  1960),  the  Gulf  coast  of 
Mexico  (Psuty,  1966),  the  Lake  Michigan  beach  (Evans,  1939),  the  Virginia  beach 
(Harrison-Wagner,  1964),  and  the  Outer  Banks  beach  (Sonu- Russell,  1966). 

TRANSITIONAL  REGIME 

This  regime  represente  the  ca6e  in  which  the  wave-current-topography 
interaction  is  undergoing  transition  caused  by  changes  in  the  wave  field.  However, 
a  systematic  description  of  the  current  Held  seems  possible  as  long  as  the  topo¬ 
graphy  remains  relatively  stable.  An  example  is  shown  in  Figure  13  (Shadrin, 
1961).  The  current  velocities  are  generally  greater  near  the  embayment  than 
near  the  shoal,  and  the  rip  position  is  associated  with  the  angle  of  wave  incidence. 
Similar  results  have  been  reported  by  Riviere  et  al  (1961)  on  the  French  Mediter¬ 
ranean  coast,  Hom-roa  and  Sonu  (1962)  on  the  Pacific  coast  of  Japan,  and  the  en¬ 
gineers  of  Nii|£da,  Japan  (Figure  14). 

As  the  obliquity  of  wave  incidence  increases,  the  circulating  cells  become 
elongated  parallel  to  the  shoreline  and  take  on  the  patterns  such  as  observed  on 
the  Outer  Banks.  As  already  described,  the  rip  pulsation  gives  way  to  a  more 
continuous  outflow,  issuing  diagonally  toward  offshore  and  being  displaced  gra¬ 
dually  downstream  with  the  migration  of  the  sand-wave  system. 

Under  natural  conditions,  the  Interaction  between  waves,  curren*  s,  and  the 
topography  is  believed  to  be  a  transition  from  one  quasi-equilibrium  to  another. 
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Consequently,  the  transitional  regime  may  represent  the  most  frequently  encoun¬ 
tered  case  in  the  field. 

FORCED  EQUILIBRIUM  REGIME 

Man-made  structures  or  fixed  topographic  features  may  impose  a  certain 
constraint  on  the  variables  associated  with  the  equilibrium  regime.  This  seems 
to  be  best  typified  by  the  case  reported  from  the  Scripps  beach,  California,  where 
topographic  diversification  in  the  vicinity  of  the  submarine  canyons  forces  inci¬ 
dent  waves  to  refract  into  a  limited  number  of  convergence  -  divergence  patterns 
(Shepard -Inman,  1950).  As  a  result,  an  extremely  steep  gradation  of  the  along¬ 
shore  breaker-height  distribution  develops  on  this  beach  -  namely  on  the  order 
of  1:2  or  even  1:3  per  less  than  half  a  mile  of  shoreline  length  -  which  will  be 
rarely  encountered  on  ordinary  beaches  having  relatively  smooth  contours  on  the 
offshore  bid.  Since  the  long-period  wave  responds  most  sensitively  to  a  forced 
refraction,  the  nearshore  water  circulation,  reported  by  Shepard  and  Inman 
(1950),  is  most  pronounced  during  the  activity  of  long-period  waves.  The  constra¬ 
int  on  the  wave  refraction  will  also  give  rise  to  an  areal  fixation  of  the  current 
pattern  as  well  as  of  the  relief  patterns  close  to  the  shore.  Munk  and  Traylor 
(1947)  have  reported  that  a  rip  current  is  frequently  observed  at  the  head  of  a 
submarine  canyon.  This  appears  to  be  one  of  the  significant  distinctions  from  the 
natural  equilibrium  regime,  in  which  the  rip  position  is  considerably  varied  aB 
the  sand  waves  migrate  In  the  direction  of  the  predominant  longshore  currents. 

The  basic  hydraulic  mechanism  governing  the  forced  equilibrium  regime  thus  is 
likely  to  be  similar  to  that  of  the  natural  equilibrium  regime,  consisting  mainly  of 
the  processes  of  momentum  transfer  over  the  bar  crest  and  mass  transfer 
along  the  shoreline. 

SUMMARY  AND  CONCLUSIONS 

The  results  of  the  analyses  are  summarized  as  follows: 

1.  Approximation  of  the  mean  longshore  current  velocities  using  the  mul¬ 
tiple  regressive  scheme  of  both  linear  and  non-linear  types  seems  to  be  substan¬ 
tially  influenced  by  the  seasonal  and  regional  characteristics  of  the  sample  popu¬ 
lations  analyzed. 

2.  Of  all  the  independent  variables  associated  with  the  mean  alongshore 
component  of  longshore  current  velocities,  the  angle  of  wave  incidence  distingui¬ 
shes  itself  as  the  most  influential  variable  in  both  the  linear  and  the  non-linear 
product- form  regressions. 

3.  The  analytical  schemes  based  on  simplified  conditions  result  in  either 
an  undefinable  scatter  (Putnam- Munk- Traylor,  1949;  Inman-Quinn,  1951;Nagai, 
1954),  or  systematic  deviations  (Brebner-Kamphuis,  1963;  Galvin- Eagleson.  1965) 
from  field  observations.  It  may  be  that  the  mechanism  of  the  phenomenon  has  not 
been  truthfully  accounted  for  in  these  schemes,  or  that  the  observations  failed  to 
produce  an  unbiased  mean  of  the  multiple-velocity  current  field.  Namely,  the 
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friction  coefficient  in  the  Putnam -Munk- Traylor  formula  may  require  further  re¬ 
finement  by  taking  into  consideration  the  dynamic  processes  of  energy  dissipation 
in  the  surf  zone  environment. 

4.  The  temporal  variability  of  the  current  velocities  appears  to  be  small 
as  compared  with  the  spatial  variability.  The  temporal  variability  of  both  the 
current  speeds  and  directions  appears  to  be  more  pronounced  over  the  offshore 
bar  than  inshore  of  this  position. 

5.  The  spatial  variability  of  longshore  currents  -  including  the  processes 
of  acceleration,  decelaration,  stagnation  reversal,  and  occasionally,  bifurcation 
-  arises  from  the  fact  that  the  current  field  is  composed  of  a  multitude  of  velocity 
vectors  whose  distribution  is  mainly  influenced  by  the  regimes  of  the  wave-current- 
topography  Interaction  in  the  surf  zone.  Thus,  the  longshore  current  may  be  clas¬ 
sified  into  four  major  types  of  velocity  field  in  accordance  with  these  interaction 
regimes  -  the  hypothetical  regime,  the  natural  equilibrium  regime,  the  transi¬ 
tional  regime  and  the  forced  equilibrium  regime. 

6.  The  current  field  of  the  natural  equilibrium  regime  appears  to  be  gene¬ 
rated  by  the  momentum  transfer  from  plunging  breakers  over  the  bar  and  the 
mass  transfer  from  spilling  breakers  over  the  shoal  as  well  as  from  swash  acti¬ 
vities  on  the  subaerial  slope.  The  seaward  mass  discharge  by  rip  currents  takes 
care  of  the  hydrostatic  potential  resulting  from  the  latter. 

7.  Under  natural  conditions,  the  nearshore  topography  participates  in  the 
longshore  current  mechanism  as  a  dynamic  variable,  not  only  redistributing  the 
breaker  influx  into  different  positions  along  the  shore  but  also  itself  undergoing 
displacements  and  transformation  due  to  the  waves  and  the  currents  thus  affected. 
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